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Abstract - We introduce a new model of hard spheres under confinement for the study of the glass 
and jamming transitions. The model is an one-dimensional chain of the d-dimensional boxes each 
of which contains the same number of hard spheres, and the particles in the boxes of the ends 
of the chain are quenched at their equilibrium positions. We focus on the infinite dimensional 
limit (d ^ oo) of the model and analytically compute the glass transition densities using the 
replica liquid theory. From the chain length dependence of the transition densities, we extract the 
characteristic length scales at the glass transition. The divergence of the lengths are characterized 
by the two exponents, —1/4 for the dynamical transition and —1 for the ideal glass transition, 
which are consistent with those of the p-spin mean-field spin glass model. We also show that the 
model is useful for the study of the growing length scale at the jamming transition. 


Introduction. — When liquids are rapidly cooled well 
below the melting temperature, the relaxation time and 
viscosity drastically increase and eventually they freeze 
into the disordered state. This phenomenon is known as 
the glass transition [1]. The mechanism of this slowing 
down and the presence/absence of an underlying genuine 
phase transition are under active debate Hi]. 

A good way to understand theoretically the glass tran¬ 
sition would be to first construct a mean field theory and 
then consider finite dimensional effects. Recently the glass 
transition of the infinite dimensional hard spheres is stud¬ 
ied as a mean-field model, and the exact thermodynamic 
theory of the model is developing [SHE]- The glass transi¬ 
tion of the model can be characterized by the two relevant 
densities [5]. The first is the dynamical transition den¬ 
sity, <pd, below which the phase space is divided into an 
exponentially large number of metastable states. In the 
mean-field limit, the dynamics of the system gets frozen 
at this density as the system is trapped into one of the 
metastable states. With increasing the density, the num¬ 
ber of the metastable states decreases and eventually be¬ 
comes sub-exponential, at which the system undergoes the 
thermodynamic transition into the ideal glass state. This 
transition density tpK is called the Kauzman density. 

What types of finite dimensional effects play a dom¬ 
inant role in the glass transition is not clear yet mm- 


One promising theoretical scenario which incorporates an 
effect of fluctuations of finite dimensional systems is pro¬ 
posed as the Random First Order Transition (RFOT) the¬ 
ory m- This theory assumes that the divergence of the 
relaxation time at g?d is avoided, and at ip between pd and 
Pk, the liquid state can be seen as a patchwork of the local 
metastable configurations m- The theory predicts that 
the characteristic size of these domains increases with in¬ 
creasing the density and eventually diverges at pK, where 
the thermodynamic glass transition occurs. 

In practice, this length can be measured through 
the thermodynamic behavior of liquids under confine¬ 
ment m- We prepare an equilibrium configuration of 
particles and pin only the particles outside of a cavity of 
the size L. One can define and compute the critical cavity 
size at which the boundary affects the thermodynamics of 
the system in the cavity, which is referred to as the point- 
to-set (PS) correlation length. Numerical simulations of 
realistic liquid models under this confinement have been 
carried out [IMl], which confirmed that the PS length 
indeed grows as the systems slow down. However, com¬ 
puting theoretically the PS length of realistic liquids is 
still very challenging [13 [S]. One of the difficulties is 
due to the inhomogeneous nature of the confined liquids, 
which requires the inhomogeneous version of the liquid 
state theory m The p-spin mean-field spin glass model 
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under confinement has been extensively studied [151 [T^ . 
The studies showed that the PS length of the model di¬ 
verges at the Kauzuman transition point with the power- 
law behavior with the exponent —1 and another length 
characterizing the local stability of the metastable states 
diverges at the dynamical transition point with the expo¬ 
nent — 1/4 [18] . 

The main purpose of the present letter is to propose 
a new model of hard spheres under confinement, which 
can be analyzed using the homogeneous liquid state the¬ 
ory |2Q]. Focusing on the large dimensional limit of the 
model, we compute analytically the PS length using the 
replica liquid theory. 

Since the model proposed is not a spin model but hard 
spheres, the model undergoes another type of the phase 
transition, called the jamming transition, at the higher 
density This is the transition of the liquid state 

into the state where hard spheres are packed so closely 
that the pressure diverges. We additionally show that the 
model can be used to study the growing lengths at the 
jamming transition. 

Model. — Consider an one-dimensional chain of L -I- 2 
boxes. Each box is a d-dimensional cube of the volume 
V and contains N hard spheres. We assume that a parti¬ 
cle interacts with the other particles in the same box and 
those in the nearest two boxes. Therefore, the Hamilto¬ 
nian of the model is 

L+l N L N 

1=0 i<j 1=0 ij 

where the d-dimensional vector denotes the position of 
the i-th particle in the l-th box and v(r) is the interaction 
potential between hard spheres of the diameter D, hence 
v(r) = oo (r < D), 0 (r > D). We study the model under 
the condition that the particles in the 0-th and [L + l)-th 
boxes are quenched at their equilibrium positions. Thus 
the model is interpreted as the d-dimensional hard spheres 
confined by the amorphous walls. We sketch a typical 
configuration of our model schematically in Fig. [T] 

As it is clear from the Hamiltonian, eq. ([TJ , we can treat 
the model as a L -I- 2 component homogeneous liquid, which 
enables us to analyze the model using the homogeneous 
liquid state theory. This is one of the advantages of this 
model. In particular, the free energy of this model can 
be evaluated analytically at the large dimensional limit 
using the replica liquid theory. In the following sections, 
we calculate the two relevant densities ipd and (fK for the 
glass transition as a function of the number of the boxes 
L. By interpreting L as the correlation length, we extract 
the lengths and which diverge at (fd and tpx, re¬ 
spectively. Additionally, we also apply this method to the 
jamming transition to discuss the possibility to compute 
its growing length scales. 



Fig. 1; A typical configuration of the model: The position of 
the particles in each boxes labeled 0-th to L-|- l-th are schemat¬ 
ically shown as cycles. The shaded cycles represent the frozen 
particles in 0-th and {L-\- l)-th boxes, and the filled cycles rep¬ 
resent the mobile particles in each box. In the l-th box, the 
particle positions in the (Z — l)-th and (Z -|- l)-th box are also 
shown as the solid and dashed cycles respectively to indicate 
the particles in Z-th box can not overlap with the particles in 
(Z — l)-th and (Z -|- l)-th boxes. 


replica liquid theory assuming the one-step replica sym¬ 
metry breaking (IRSB) ansatz [5l[22l[23]. The main idea 
of the theory is to consider m copies (replicas) of the orig¬ 
inal system. There is no interaction between the replicas, 
and therefore the partition function of the replicated sys¬ 
tem is written as 

m 

Z^ = \[ (2) 

a—1 


where H is the Hamiltonian of the original system, and T 
is the temperature. From now on, we set T = 1 because 
the temperature is an irrelevant variable for hard spheres. 
Within the replica theory, the logarithm of the number of 
the minima, which is called the configurational entropy, 
E, can be calculated by the formula [23] : 


/ log \ 
N dm \ m J 


(3) 


If one takes the large dimensional limit (cZ —>■ oo) and 
the thermodynamic limit (iV —>■ oo, E —>• oo, AZ/E=const), 
the calculation of the partition function of the replicated 
system can be greatly simplified because the higher order 
terms of the Mayer cluster expansion become negligible, 
and only the first term needs to be considered [51150] . The 
replicated free-energy of our model is 

log Zm = X! / (1 “ log Pi(®)) 

1=1 

+ 5 ^ / d^(Wpi{x)pi{y)f{x-y) 

^ 1=1 •' 

L — 1 „ 

+ X] / ^(^Pi(x)pi+iiy)fix -y) 

1=1 •' 

+ J dxdypi{x)po{y)f{x - y) 


Free energy at large dimensional limit. — To 

compute the free energy of the model, we employ the 


dxdypL{x)pL+i{y)f{x - y), 


(4) 
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where x = denotes the positions of 

hard spheres in the replicated space, and pi{x) = 

(j2i Y\^=i ~ is the density distribution of the 

1-th box. The boundary conditions are included in the 
last two lines of eq. where po{x) = S{x — x^) and 
Pl+i{x) = ~ microscopic density 

distributions of the frozen particles belonging in the 0 -th 
and {L + l)-th boxes. and xf~^^ are taken from their 
equilibrium positions. If the system is sufficiently large, 
the self-averaging properties for the free-energy holds: 

log Zjn ^ log Zm, (5) 


where the over line denotes the average of po{x) and 
PL+iix). Since the free-energy eq. (|3]) depends linearly 
on those variables, we have only to replace as po{x) —>■ p 
and pl+i(x) —>• p, where p = N/V is the number density 
of each boxes. 

In the next step, we introduce the IRSB Gaussian 
ansatz [^: 

/ m 

dXW_"fAi{Xa-X), ( 6 ) 

a=l 


where 7 ^( 2 :) = exp (—a;^/27l)/(27rA)'^/^. This ansatz 
claims that the distribution of the replicated par¬ 
ticles is the Gaussian with the variance Ai = 
/m(m — 1). Substituting eq. (jH) into 
eq. ® and taking the large dimensional limit [5], one ob¬ 
tains the analytical expression of the free energy: 


logZrr 

N 


— ^id T ^int 
L 


s.d = Y. 


Sint — 


3Lip 


. . Ai . m 
(m - 1 ) log -TT + log — + m 
d 




2 


1=1 


+ ‘P 


Ai + Ai+i 


1=0 


(7) 


where we introduced the normalized volume fraction ip = 
2‘^ip/d and the normalized cage size Ai = Aid^/D'^. 
GmiA) is the auxiliary function given by 


GmiA) 



0 


( 

V \/ 4 l 


m 


^(y) 


( 8 ) 


where 0(a;) = [erf(a:) -|- 1] /2. In eq. (O, the cage size of 
the 0-th and {L + l)-th boxes are zero, Aq = Al+i = 
0 , because of the boundary condition. Ai for other I is 
calculated from the saddle point conditions for Ap. 


J_ ^ 2 dS,nt 

Ai (1 - rn)d dAi 


(1 = 1,...,L). (9) 


Substituting the solution of this set of equations into 
eq. 0 , the free energy of the m replicated system is ob¬ 
tained. 



Fig. 2: The density dependence of the cage size for several L’s. 
The solid line is for L = 2, the dashed line is for L = 4 and the 
dotted line is for L — 8. 


Glass transition. — Here, we evaluate the correlation 
length of the glass transition. As mentioned before, there 
are the two relevant densities for the glass transition, ipd 
and pk- Their values can be computed by analyzing the 
free energy at m = 1 . 

First, we examine the behavior of our model near the 
dynamical transition point, pd: where the exponentially 
many metastable states emerge on the free energy. In the 
framework of the replica liquid theory, the order parameter 
to characterize the metastable states is the cage size, which 
is given by the solutions of eq. ([9]) [5]. Their numerical 
solution for Ai of the T/ 2 -th box, which is the farthest 
from the boundaries, are shown in Fig. [51 For the small 
densities, the cage size Ai ^/2 is infinity Ai^/ 2 ) = 0 ); 

meaning that the system is ergodic and in the liquid phase. 
Increasing the density, the cage size jumps discontinuously 
to a finite value at the dynamic transition point, pd- 

From Fig. [51 it is clear that pd increases as increasing 
L. From this result, we can convert Pd{L) to L{(p), the 
characteristic number of the boxes as a function of the 
transition density. We define the correlation length by 
^d = dj[(p)l2 and plot it in Fig. [3] [18]. We find that 
diverges at pd as ^d ~ {pd - p)~^l^. 

Next, we evaluate the correlation length near the Kauz- 
man density by calculating the configurational entropy, S, 
which is obtained by plugging the free energy into eq. ©• 
The final expression in the large dimensional limit becomes 


Z,{m,p,L) 




+ 0{d), 


( 10 ) 


up to the order of O(dlogd). From this expression, it is 
clear that S vanishes at 


Pk{L) 


dlogd 

3(l-b2/3T)’ 


( 11 ) 


Following the same argument as for ^d, the growing length 
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Fig. 3: The correlation legnths near the tpd and ipth- The 
filled cycles are for as a function of e = {(pd — p)lpd- The 
filled squares are for ^th as a function of e = {pth — p)lpth- 


around the Kauzman density is 


=iccp 


L 

"2 


- (p) ’ 


( 12 ) 


where px = dlogd/S is the Kauzman density of the bulk 
system (L —oo). 

Jammintg transition. — Hard spheres also undergo 
the jamming transition when the system is compressed 
quickly or making the pressure infinity [21] . The infinite 
dimensional hard spheres serve as a mean-field model of 
the jamming transition as well and are studied extensively 
using the replica liquid theory. The theory showed that 
the jamming transition is also characterized by the two rel¬ 
evant densities The jamming state obtained by com¬ 
pressing the ideal glass state up to the infinite pressure is 
referred to as the glass close packing, and its density is 
denoted as (pccp- If the jammed stated is prepared by a 
fast compression of low density hard spheres, the pressure 
becomes infinity at a much lower density than pccp [24] . 
The lowest density of the jammed states pth is called the 
threshold density. 

In this section, we show that the growing length scales at 
the jamming transition can be computed from the analysis 
of the jamming transition of the present model. We essen¬ 
tially follow the strategy developed in the previous section. 
Using the replica liquid theory with the IRSB ansatz, we 
compute Pth and pccp of the model as a function of the 
number of the boxes, L, which naturally gives rise to the 
two characteristic lengths ^th{p) and iacpi.p)- Note that 
recent studies showed that the full RSB ansatz is needed 
for the fully exact computation near the jamming transi¬ 
tion mi- However such a computation for the present 
model seems quite involved. Here, we wish to stick to 
the IRSB ansatz and demonstrate that the replica theory 
analysis of the model provides the analytical expressions 
of the growing length scales at the jamming transition. 


In the replica liquid theory, the divergence of the pres¬ 
sure corresponds to take the m —>■ 0 limit, because m is 
inversely proportional to the pressure [5] . Thus pth can be 
calculated by taking the m —0 limit in the self-consistent 
equation for the order parameter eq. Since the cage 
size, Al, vanishes at the m —)■ 0 limit, we set Ai = mai 
0. This is substituted to eq. (0) before the m 0 limit 
is taken. We numerically solve the equation and find that 
the behavior of q;l /2 is qualitatively the same as that of 
^L /2 for the glass transition (Fig. 0). At low densities, 
Q^l /2 = oo is the only solution of eq. (0). This means that 
there are no jammed states at those densities. At higher 
densities, the solution of eq. ® becomes finite. This tran¬ 
sition density is the lowest density, pth, of the jammed 
states. As in the case of the glass transition, the charac¬ 
teristic length Cth{p) is computed from the L dependence 
of Pth- This length ^th is also plotted in Fig. 0 We find 
the power-law divergence ^th ~ {pth — ■ Likewise, 

we analyze the correlation length near pocp- The loga¬ 
rithm of the number of the jammed states is calculated by 
the m —^ 0 limit of the configurational entropy eq. (nni 
and PGCP is defined as the density where the configura¬ 
tional entropy becomes zero. From eq. m, it is clear that 
the behavior of pccp should be the same as those of px 
because eq. m is independent from m. Thus we obtain 
the critical behavior of the diverging length near pocp as 
Cgcp ~ {pGCP - p)~^- 

Summary and discussion. — In this letter, we con¬ 
sidered an one-dimensional chain of the d-dimensional 
boxes each of which contains N hard spheres as a model 
system of liquids under confinement. By focusing on 
the large dimensional limit, we analytically computed the 
phase diagram of the model. From the chain length de¬ 
pendence of the transition densities, we derived the criti¬ 
cal behavior of the two relevant length scales of the glass 
transition. We also showed that there are the two relevant 
length scales of the jamming transition, whose critical be¬ 
haviors are similar to those of the glass transition as long 
as the IRSB ansatz is assumed. One of the advantage of 
the model is that, although the system is under confine¬ 
ment, it can be regarded as a homogeneous liquid, and 
thus the model can be analyzed using the usual replica 
liquid theory. 

We summarize the results in Fig. |4l The upper panel 
is the phase diagram of the glass transition of the con¬ 
fined liquid, where p is the normalized density p = 2'^p/d, 
and L is the distance between the boundaries. There are 
the two relevant length scales and which diverge 
as ~ {pd - py^^'^ and ^x ~ {pk - , respec¬ 

tively. When L > ^d, the system holds the ergodicity 
and is in the liquid phase. For ^x < L < ^d, the phase 
space splits into the many sub-spaces and the ergodicity 
of the system is broken at least in the large dimensional 
limit. In L < ^x region, the number of the sub-spaces 
becomes sub-exponential and no longer contributes to the 
entropy. Thus the equilibrium phase transition from the 
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^th ^ {‘fth — ‘p) ^GCP ~ {‘PGCP - 0) ^ 





Fig. 4: The schematic phase diagrams of our model: The 
upper and lower panel is the phase diagram for the glass and 
jamming transition, respectively, (p = 2'^ipjd is the normalized 
volume fraction and L is the confinement length. The insets 
are schematic pictures of the phase space in each phase, where 
the cycles represent the accessible regions in the phase space. 


liquid phase to the ideal glass phase occurs. These re¬ 
sults are consistent with the results for the p-spin mean- 
field spin glass model [T8l[T9] , where the critical behaviors 
^ {T — and ^ {T — Tk)~^ are observed. 

This confirms that the proposed model is considered to be 
a hard sphere extension of the p-spin mean-field spin glass 
model under confinement. The result for is also con¬ 
sistent with the prediction of the inhomogeneous mode¬ 
coupling theory for the dynamic correlation length diverg¬ 
ing at the mode-coupling transition point HZ]. 

The lower panel of Fig. 2] is the phase diagram of the 
jamming transition. There are the two relevant length 
scales S,th and ^gcp which diverge as ^th « {0th - 
and ^GCP ~ {0GCP - respectively. When L > ^th, 

there are no jammed states. For ^gcp < L < ^th, there 
exist exponentially many jammed states. The number of 
the jammed states decreases with L, and becomes sub¬ 
exponential when L < ^gcp- 

We should emphasize that the computation of ^th and 
^GCP in this work is not exact, since the IRSB ansatz is 
recently found to give unstable solution near the m —>■ 0 
limit [7]. If we use a better ansatz, the critical exponents 
may change [8]. Regarding this point, we wish to indicate 
the two points: (1) Even though the solution is unsta¬ 
ble, the physics of the IRSB solution can still survive in 
certain situations |25j and thus the analysis given here is 


still useful 0. (2) In case of ^ggp, the results would not 
be effected by this instability at least in the large dimen¬ 
sional limit, because the expression of the configurational 
entropy is m independent. 

Finally we discuss the physical meaning of ^th and ^gcp- 
Several different types of lengths are known to diverge 
at the jamming transition and the relations between the 
lengths are still not clear enough [26H36] . Since ^th and 
Cgcp represent the distance between boundaries at which 
the jammed configurations first appear/disappear, we ex¬ 
pect that these lengths correspond to the length detected 
in the finite size scaling analysis \<p — [5^ . 

However this is a tentative expectation, because we ap¬ 
plied a specific boundary condition (amorphous bound¬ 
ary) which is suitable to compute the PS length at the 
glass transition. It is not clear which length should be de¬ 
tected by this condition for the jamming transition. To re¬ 
solve this question, it should be useful to study the present 
model under various different types of confinements such 
as the flat walls m or random boundary |19j . 
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